Dissipative Stern-Gerlach recombination experiment 
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The possibility of obtaining the initial pure state in a usual Stern-Gerlach experiment through 
the recombination of the two emerging beams is investigated. We have extended the previous work 
of Englert, Schwinger and Scully Q including the fluctuations of the magnetic field generated by 
a properly chosen magnet. As a result we obtained an attenuation factor to the possible revival of 
coherence when the beams are perfectly recombined . When the source of the magnetic field is a 
SQUID (superconducting quantum interference device) the attenuation factor can be controlled by 
external circuits and the spin decoherence directly measured. For the proposed SQUID with dimen- 
sions in the scale of microns the attenuation factor has been shown unimportant when compared 
with the interaction time of the spin with the magnet. 



I. INTRODUCTION 

The Stern-Gerlach experiment gives an experimental 
evidence of the quantum nature of the spin of a particle. 
When a beam of spin 1/2 particles, in the eigenstate |+) 
of S x , goes through a variable magnetic field in the z 
direction and the outcoming particles are detected on 
a screen, we observe the presence of two distinct peaks 
corresponding to the spins in the positive and negative z 
direction. 

Moreover, besides giving an evidence of the spin quan- 
tization, this example is considered the paradigm of the 
measurement process. The beam goes into the magnet 
in a pure state, for example, the eigenstate \+) x of S x , 
and as far as the measurement of S z is concerned it is de- 
scribed as a superposition of the eigenstates of S z . Later, 
as a result of the measurement process, the quantum 
state collapses into one of the eigenstates of the latter. 

A question that arises is: could we get the initial pure 
state of the system if we recombine the two beams from 
the Stern-Gerlach experiment? 

This idea of the recombination of the beams, through 
the so-called Stern-Gerlach Interferometer(SGI), is an old 
one, but has always been treated in a qualitative way 
One of the first quantitative studies was presented in a 
series of three articles 0, IE Q| , in which the authors ex- 
amined the possibility of recombination and concluded 
that the precision of control of the magnetic field is of 
fundamental importance to the reconstruction of the ini- 
tial pure state. However, in those articles the authors did 
not take into account the origin of the magnetic field and 
its possible sources of fluctuations. 

In a more recent study Q this problem has been revis- 
ited and the quantum nature of the magnetic field taken 
into account. However, only the coupling of the uniform 
part of field to the spin of the particle has been consid- 
ered and the effects of its fluctuations on the spatial part 
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of the wave function have been disregarded. 

In another manuscript @| a model for a dissipative 
Stern-Gerlach experiment was proposed without explicit 
reference to the origin of the fluctuations. It basically 
dealt with a Stern-Gerlach apparatus inside a dissipative 
medium and, besides, did not consider the recombination 
process. 

In this manuscript we intend to make a fully quantum 
description of the magnetic field and study the effects 
of its unavoidable fluctuations in the spatial part of the 
wave function. We also pay special attention to the origin 
of these fluctuations by presenting a model where the 
magnetic field is generated by a pair of SQUIDs where the 
magnetic flux is known to obey a Langevin equation. We 
will be only interested in the spin coherence that can be 
measured through the mean value of S x and arises from 
the off-diagonal terms of the reduced density operator of 
the system in the spin space. 

Another point that is worth mentioning here is the fact 
that we will be aiming at a problem quite different from 
the standard one in the area of dissipative systems. Usu- 
ally people want to get rid of decoherence for macroscopic 
quantum mechanical variables whereas we will try to in- 
crease the effect of decoherence on a microscopic variable. 
The reason for this is to fully test the existing theory of 
decoherence as unequivocally due to the sources of noise 
and/or dissipation known to be coupled to the system. 

The paper is organized as follows: in Sec. II we present 
a model for the SGI where the magnetic field is generated 
by a pair of SQUIDs. In Sec. Ill we evaluate the dynam- 
ics of the reduced density operator in the spin space using 
the Feynman- Vernon path integral approach. In Sec. IV 
we analyze the coherence revival in the SGI. Finally, we 
summarize our results in Sec V. 



II. THE STERN-GERLACH 
INTERFEROMETER 

The SGI can be divided in two parts: the first one is 
responsible for separating the beams and the second one 
for recombining them. We will follow the basic model and 
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Figure 1: Loss and revival of coherence in a SGI without 
fluctuations. 



approximations of Q, considering that the beams split in 
the z direction due to the magnetic field gradient which 
also points along this direction. As referred to a given 
origin, half way from two specific magnets (see arrange- 
ment below), the magnetic field should depend on y in 
an anti-symmetric way to enable a perfect recombination 
of the two beams. In order to solve the problem we will 
consider that the velocity in the y direction is constant 
which allows us to replace its position dependence by a 
time dependence. 

Our model Hamiltonian is 



H = ^—-a z f(t)z, 
Am 



(1) 



where f (t) = fidB/dz and /i is the magnetic moment 
of the particle. Here we are not considering the uniform 
part of the magnetic field that is only responsible for the 
precession the spin about the z direction and has already 
been treated in Q. 



For a perfectly noiseless anti-symmetric field with a 
piecewise constant gradient we observe loss and revival 
of the spin coherence through the behavior of the off- 
diagonal elements of the reduced density operator in the 
spin space, as shown in Figure 1. There, T is the exper- 
iment time and the revival in the middle of the experi- 
ment is due to the recombination of the wave packets in 
the momentum space. The vertical lines show four dif- 
ferent regions of a SGI. In the first part the beams are 
accelerated in opposite directions along z. The second 
and third parts are due to the empty interval between 
the magnets where the direction of the acceleration for 
each beam suddenly changes. In the last part the acceler- 
ation of each beam returns to its original value and they 
are finally recombined as they emerge from the set up. 



We now introduce the fluctuations of the magnetic field 
to the problem. In order to do this we will consider a 
SGI where each of the two magnets is formed by a pair 
SQUIDs arranged as shown in figure 2 . With this model 
we obtain a field configuration very similar to that pro- 
posed in 2| and described above. 





Figure 2: Lateral view of the proposed apparatus for the ISG and the magnetic field lines in a frontal view. 



The total flux inside a SQUID has its dynamics de- 
scribed by the model Hamiltonian given by Q 
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function J (lu) = 

2 7^UZ U ~ ^k) - rju. Here $' is the fluctu- 
ation about one of the many metastable flux values 



at zero external field. We are considering that there 
is no tunnelling or thermal activation of the flux 
variable to its neighboring minima. Thus, the flux is 
oscillating around a local minimum close to n$o and 
Lq = [1/L + 2irio/n$o] is an effective inductance, 
that arises from a second order expansion around the 
minimum of the electromagnetic potential energy, zo is 
the critical current for the SQUID. For this approxima- 
tions to be valid we should be sure that we have many 
minima, in other words: 27rioi/$o 3> 1 ■ As a matter 
of fact we shall assume that each SQUID is carrying 
a persistent current corresponding to the metastable 
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minimum at n = 1. 

Here, it should be emphasized that we are only taking 
into account the parameters of one single SQUID of each 
pair, say, the smaller one. Operationally it means that 
the particle trajectories are always close to that same 
circuit over half of the experiment time. In the second 
half of the experiment the same applies to the second 
pair of SQUIDs. The presence of the larger SQUID in 
each pair is only to mimic the field lines of a usual Stern- 
Gerlach apparatus. We believe that this requirement can 
be dropped at the expense of using effective parameters 
for the coupled SQUIDs in a Hamiltonian of the same 
form as (0). 

Therefore, we are dealing with a particle coupled to a 
magnetic field that is produced by a flux which, on its 
turn, is coupled to a bath of harmonic oscillators. Now 
if we write the magnetic field in terms of the flux we will 
have the particle indirectly coupled to the bath through 
the flux. This can be done if we write $ = B (z) A (z) 
where B (z) is the magnetic field at z as seen by the 
particles whereas A (z) is an effective area crossed by 
the field lines through which the magnetic flux is exactly 
given by the total value of <£> inside the ring. Since the 
latter does not depend on z we have 



, . dA . . . dB 
B{z) — + A{z) — 
oz oz 

that can be manipulated to give 



(3) 



dB 

dz 



a (z) $ , with a (z) = 



1 dA 
A{zf~dz~ 



(4) 



Using of the last expression, the total Hamiltonian 
reads 
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For an original Ohmic spectral function, which is 
known to hold for SQUIDs, we obtain the following ef- 
fective spectral function: 



with e = /ia . 

Now we have the particle coupled indirectly to the bath 
of oscillators. Following the prescription of || we can 
eliminate the flux variable in the Hamiltonian and write 
a new one where the particle is coupled directly to new 
oscillators with an effective spectral function. The new 
Hamiltonian is 
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with r) = e 2 Ll/R, ft = l/s/T%/W - 2CL and ft' = 
1/^/CLq, Here R and C are, respectively, the resistance 
and capacitance of the junction of the SQUID. The ef- 
fective spectral function is Ohmic with a new cutoff fre- 
quency given by the minimum of ft and ft' . 

For typical SQUIDs we have the following values: 
C ~ 10- 12 F, L - 1CT 10 #, i - 1CT 5 A and R - in. To 
obtain a (z) we will suppose a SQUID with dimensions 
of (lCT 5 x 10~ 3 ) m 2 and will estimate the total flux us- 
ing the field in the middle of two infinite wires which is 
multiplied by the area of 10~ 8 m 2 . Comparing this with 
the flux at a distance z we can obtain an expression for 
A(z) and then estimate a(z). For z = 10~ 3 m we have 



a{z) = 10 



13„ 



These specifications for the apparatus 



give L ~ 10- 10 ff and ij = lO-^Kgs- 1 . For cooper 
atoms, m = 1,8 x 10 _25 fc<7, the relaxation time, 7 _1 , is 
of the order of 10 15 seconds. 

We end this section with this proposal for a model for 
the SGI that incorporates unavoidable sources of fluctu- 
ations of the magnetic field. In our model these fluctua- 
tions are modelled by a coupling to a bath of oscillators 
following Q. Now we have to evaluate the reduced den- 
sity operator of the system in the spin space to see the 
behavior of the spin coherence through its off-diagonal 
elements. 



III. THE REDUCED DENSITY OPERATOR 

Our model of a pair of SQUIDs generating the mag- 
netic field for the SGI obeys the following Hamiltonian 



H — Hq + Hj + Hf, 



with 



(5) 
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(9) 
(10) 

.(11) 



It describes the dynamics of a particle submitted to a 
linear potential and also coupled to a bath of oscillators 
both in a spin dependent way. The problem of a particle 
in a linear potential and coupled to a bath of oscillators 
has already been treated in the literature and the only 
difference here is the spin dependence that was absent in 
@. However, since we have only one spin operator in JHJ 
it acts as a parameter in our problem and we will use a 
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slightly modified Feynman- Vernon formalism to solve it. space are written as 
We are not going to show all the details, of the calculation 
for which the reader should consult fft. IflL HH] . 

For separable initial conditions, where the inter- 
action between the system and the bath is turned , , _ f f , , , f , , ( , , f , 
on at i = 0, one has (a/R'S] p (0) \y'Q'S') = Pss' [x,y,t) - J J dx dy J ss , [x,y,t,x ,y , U) {x b\ p x (Uj \y b ) 

(x'S\ p x (0) \y'S') (R'| P 2 (0) |Q') and the four elements of (12) 
the reduced density operator of the system in the spin with 



J ss , [x, y, t; x', y',0)=J J J dRdR'dQ'K a (x, R, t; x', R, 0) (R'| p 2 (0) |Q') K* s , (y, R, S', t; y', Q', 0) (13) 

I 

which propagates the reduced density operator in time. (N-dimensional vectors) of the bath of oscillators. Within 
In this expressions R, R', Q' are arbitrary configurations the path integral formalism it is written as 

I 



J ss ,(x,y,t;x>,y',0)= f Dx(t') f Dy (t>) e *( s o'KOMS'[»C)]) Fss , [ x tf) , y (f)] ( 14 ) 

J x' J y' 



with 



F ss ,[x(t'),y(t')} = 



J J J dRdR'dQ'p 2 {R',Cl',0) J DR{t') J Q Dq,{t') 

xe i(s|[ a5 (t'),R(t')]-S!'[i/(t').Q(t')]) e *( s R[ R ( t ')]- s B [«(*')]) (15 ) 
I 



being the so-called influence functional, which contains in which 
all the influence of the bath on the system. It has been 
evaluated before and the resulting expression is 

F ss ,[x(t'),y(t')] = e-^*A*nv(t')}, (16) 



$ ss ,[x(t'),y(t')} = —[ Sx (0)-S'y(0)} I dt' 7 (t')[Sx (t') - S'y (t')} + 



Jo 

+ ^ / dt ' f dt " t 5x (*') - S 'y (*')] 7 (*' - t") [Si (t") + S'y (t")] + 

+ f dt' [ dt" [Sx (f) ~ S'y (t')} a R (t' - [Sx (t") - S'y (t")} , (17) 
Jo Jo 



r 

with 



j(t) = —9(t) [ ^1 cos cut duj (18) 

TO7T J LU 
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and 



we have an Ohmic case, eq [7], given by 
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a R (f - t") = -J J(u>) coth \—^, ) cos u {f - t") du: 



T(, ,\ - / ^ f° rUJ < ^ 
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(20) 



(19) 



All we have to know is the spectral function of the where we have introduced a cutoff frequency ft. With 



bath, 



this we obtain 
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J ss ,(x,y,t;x',y',0) = f Dx(t') P ' Dy{t') e^'K'X*' 

J x' J y' 

(21) 



that completely characterizes it. In our proposed model with 
I 



p ss , [x (f) , y (*')] = i J dt' (x 2 - y 2 ) - f (t') (Sx - S'y) ~\{Sx- S'y) (Sx + S'y) 

- f dt' f dt" [Sx (t') - S'y (t')] a R (t' - t") [Sx (t") - S'y (t")} . 
Jo Jo 

Computing the diagonal (S = S' = ±1) and the off-diagonals (S ^ S' = ±1) terms we obtain 
fa [x (f) , y (t')} = i J* dt' (x 2 y 2 ) T fo (f) (x - y) - \ (x - y) (x + y) 



+ 



dt' / dt" [x {t') - y (t')\ a R (f - t") [x (t") - y {t")\ 
o Jo 



and 



(i od [x it') , y (t')\ = i J* dt' [| (x 2 - y 2 ) T fo (f) (x + y)-\{x + y) (x y) 

- { dt' f dt" [x (f) + y (t')} a R (f - t") [x (t") + y (t")} . 
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So we have a different dynamics for the diagonal and off- 
diagonal terms. Defining new coordinates q = (x + y) /2 
and £ = x — y we decouple the x (t) and y (t) trajectories 
and get 



■od = I dt' 
Jo 



mqi - rj£q T 2/o (*') q 



and 



(22) 
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t t' 

J«o*>to,t,t;M,0) = f q Dq(t') /W) (25) = dt 'j Q dt"q(t')a R (t'-t")q(t"). (29) 

J a' Jf 



with 



Evaluating the integral in the usual way we expand 
it about the classical trajectory and solve the remaining 
functional integral to obtain 



P' d = / dt' mqi-rtqTfoit')^ 



, , 1 L' a' class r>ll class] , , 

Jd(od)(q^,t;q',t',0) = e^ t3 ^ ""<«<««> J G (q, £, t; <?',£', 0) 

(30) 



(26) 



>d(od) 

with 



0d = [ dt' [ dt"Z(t')a R (t'-t")Z(t"), (27) 
Jo Jo 
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2-kH sinh 
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^ class = Z qL _( t )+ML + (t)-tfN(t) 

-?qM(t)TX(t)ZTZ(t)g, (32) 



L± (t) = m7 [coth (7*) ± 1] 



(39) 



#d - (*) + C'q'L+ (t) - (t) 

-Z'qN(t)T2X(t)qT2Z(t)q' , (33) 



TV (t) = 7717- 
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sinh7t ' 
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and 



'// class 



= -[eA(t) + 2tf'B(t)+Z' 2 C(t)], (34) 
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(41) 



class = 2 [g 2 A (t) + 2qq , B (t) + ? /2 c (t) ] 



(35) 



In the above expressions we have used £ (0) = 
£ (*) = £ an d the same notation holds for q. We have 
also defined the following functions 



A(t) 



-2ft 
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X(t) = 
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and 
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(43) 



For an initial Gaussian packet centered at the origin, 



ip(x) = 



1 



e i° 2 , 



(44) 



we have 



- 7 t 



B(t) = - 2 / dt' / ^"e 7 (*' +t ")sinh7t' 
sinh jt Jo Jo 

xa H {t 1 -t")smh-i{t-t") , 



(45) 



The diagonal element of the reduced density operator 
(38) for the system becomes 



Pd (q,£,t) 
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whereas the off-diagonal element is 
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So we have obtained an analytical expression for the re- 
duced density operator of the system in the spin space, 
when the initial spatial part of the wave function is a 
Gaussian centered at the origin. Now we are going to 
analyze this result. 



IV. SPIN COHERENCE IN THE SGI 

The positions of the center of the packets are given by 
the probability density 



Pd (q,0,t) = 



1 



27TCT (t) 



exp - 



2a (ty 
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qT M 



(49) 



r hl 




with 



Figure 3: Behaviour of the exponent of (53), showing the 
possible recovery of coherence. T , the experiment time, is 
here of the order of 10 -9 s. 



a(t) 



hy/2a (t) 
M(t) ■ 



(50) 



So, as expected, we have a Gaussian packet with width 
a (t) whose center is described by 



z{t) 



z(t) 

M(t) 



-2j(t-t') 



(51) 

which is the classical trajectory of a particle subject to a 
force /o (t) and to the viscous force rjv. 

To see how pure the state is as it evolves in the SGI we 
will investigate the off-diagonal elements of the density 
operator in the spin space. In the high temperature limit, 
KT > fry, the functions A(t), B (t) and C (t) can be 
easily evaluated and yield 



Pod{t) = h(t)exp 



1 / AS 
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M(t) 



2H^a (t) a' (t) 



and 



(53) 



htfL 




Figure 4: Factor h (t) showing the irreversible loss of coher- 
ence. T is the same as above. 
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(54) 



In the non-dissipative limit, 7 — > 0, h (t) goes to one 
and the exponential recovers the non-dissipative expres- 
sion for the reduced density operator. Therefore, we can 
say that the irreversible loss of coherence is in the factor 
h (t) . In figures 3 and 4 we show the qualitative behavior 
of the exponential term and the h (t) factor. 

In these graphs we can see that the recombination 
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T4 

Figure 5: Behaviour of the non-diagonal element in the spin 
space. 

could lead to a revival of coherence, but h (t) continually 
decreases which makes this revival negligible as shown in 
figure 5, where, in the end, we have only 20% of coherence 
left. 

Despite the qualitative nature of the graphs, we can 
say that in order to have a significant revival of coher- 
ence the time of the experiment should be much shorter 
than the decoherence time, which can be extracted from 
the expression for h (t) . Since this expression is very cum- 
bersome, we have to do it grafically. 

For the proposed model of a magnetic field generated 
by a pair of SQUIDs with dimensions of microns, we ob- 
tained 7 _1 = 10 15 s. These values and the other parame- 
ters already used give us the decoherence time, r, of 10 5 s 
for temperatures close to 0.1K. 

In the usual Stern-Gerlach experiment the particles 
have velocities of about 1000 m/s that give us an interac- 
tion time of the order of 10~ 6 s, which is smaller than the 
decoherence time. So, for particles with velocities of the 
order of 1000 m/s decoherence can not be observed and 
will not be important in the SGI. Therefore, in order to 
control the spin coherence time in an observable way we 
should use SQUIDs of smaller dimensions which will fur- 
nish us with shorter relaxations times. Larger SQUIDs 



would only make matters even worse. 

Although we have only stressed the possible effect of 
field fluctuations on the recombination one should bear 
in mind that the perfect recombination is by itself a very 
difficult task to perform. 



V. SUMMARY 

In this work we have investigated the effect of spin 
coherence, discussing experiments with beams of parti- 
cles subject to magnetic fields. We have extended the 
analysis of the Stern-Gerlach interferometer adding un- 
avoidable fluctuations of the magnetic field through the 
coupling of the system to a bath of harmonic oscillators. 
The introduction of these fluctuations adds a factor to 
the off-diagonal elements of the reduced density operator 
in the spin space that decays in a time r, the decoherence 
time. However we have shown that even for a very "noisy" 
environment provided by particularly chosen parameters 
for a SQUID the decoherence time is still extremely long 
compared to the experiment time. 

When the magnetic field is generated by our proposed 
micrometric SQUID, we observed that dissipative effects 
are not important when the particles velocities are of the 
order of 1000 m/s. We believe that changing the SQUID 
parameters it will be possible to observe the loss of spin 
coherence by reducing the SQUID size. In this case we 
will have a model of a reservoir for a genuine microscopic 
quantum mechanical variable which, in our case, is the 
magnetic moment of an atom! 
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